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Gluon field distribution in baryons 
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Methods for revealing the distribution of gluon fields within the three-quark static-baryon potential are pre- 
sented. In particular, we outline methods for studying the sensitivity of the source on the emerging vacuum 
response for the three-quark system. At the same time, we explore the possibility of revealing gluon-field distri- 
butions in three-quark systems in QCD without the use of gauge-dependent smoothing techniques. Renderings 
of flux tubes from a preliminary high-statistics study on a 12'^ x 24 lattice are presented. 



1. INTRODUCTION 

Recently there has been renewed interest in 
studying the distribution of quark and gluon 
fields in the three-quark static-baryon potential. 
While early studies were inconclusive im- 
proved computing resources and analysis tech- 
niques now make it possible to study this system 
in a quantitative manner "2"6' . In particular, it is 
possible to directly compute the gluon flux distri- 
bution |4I5| using lattice QCD techniques similar 
to those pioneered in mesonic static-quark sys- 
tems [STTIS] . 

Like Okiharu and Woloshyn 5 our first inter- 
est is to test the static-quark source-shape depen- 
dence of the observed flux distribution, as repre- 
sented by correlations between the quark posi- 
tions and the action or topological charge density 
of the gauge fields. To this end, we choose three 
different ways of connecting the gauge link paths 
required to create gauge-invariant Wilson loops. 
In the first case, quarks are connected along a T- 
shape path, while in the second case an L-shape is 
considered. Finally symmetric link paths approx- 
imating a Y-shape such that the quarks approxi- 
mate an equilateral triangle are considered. The 
latter is particularly interesting as the probability 



of observing A-shape flux-tubes are maximized in 
this equidistant case. 

Because the signal decreases in an exponential 
fashion with the size of the loop, it is essential 
to use a method to enhance the signal. We use 
strict three-dimensional APE smearing [H] to en- 
hance the overlap of our spatial-link paths with 
the ground-state static-quark baryon potential. 
Time-oriented links remain untouched to preserve 
the correct static quark potential at all separa- 
tions. 

Excellent signal to noise is achieved via a high- 
statistics approach based on translational sym- 
metry of the four-dimensional lattice volume and 
rotational symmetries of the lattice as described 
in detail below. This approach contrasts previous 
investigations using gauge fixing followed by pro- 
jection to smooth the links and resolve a signal in 
the flux distribution |4I5| . 

2. WILSON LOOPS 

To study the flux distribution in baryons on the 
lattice, one begins with the standard approach of 
connecting static quark propagators by spatial- 
link paths in a gauge invariant manner. APE- 
smeared spatial-link paths propagate the quarks 
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from a common origin to their spatial positions. 

The smearing procedure replaces a spatial link, 
Ufi{x), with a sum of 1 — a times the original 
link plus a/4 times its four spatially oriented sta- 
ples, followed by projection back to SU{3). We 
select the unitary matrix U^^ which maximizes 
TZetr{U^^ U'^}), where U'^ is the smeared link, by 
iterating over the three diagonal SU{2) subgroups 
of SU{3) repeatedly. We repeat the combined 
procedure of smearing and projection 10 times, 
with a = 0.7. 

Untouched links in the time direction propa- 
gate the spatially separated quarks through Eu- 
clidean time. For sufficient time evolution the 
ground state is isolated. Finally smeared-link 
spatial paths propagate the quarks back to the 
common spatial origin. 

The three-quark Wilson loop is defined as: 
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(1) 



where Uj is a staple made of path-ordered link 
variables 



Ql 



Figure 1. Projection of the T-shape path on the 
xy-plane. 



Figure 2. Projection of the L-shape path on the 
xy-plane. 

Ql 




Uj = Pexp ig / dx^Af^^x) 



(2) 



and Sj is the path along a given staple. 

In this study, we consider two-dimensional 
spatial-link paths. Figures 1 and 2 give the pro- 
jection of T and L shape paths in the xy plane. 
Figure 3 shows the idea behind the Y-shape. In 
all cases the three quarks are created at the ori- 
gin, O (white bubble), then are propagated to the 
positions Ql, Q2 or Q3 (black circle) before be- 
ing propagated through time and finally back to 
a sink at the same spatial location as the source 
(O). 

For the Y-shape we create elementary diago- 
nal "links" in the form of boxes as shown in Fig. 
4. The 1x1 and 1x2 boxes are the average 
of the two path-ordered link variables going from 
one corner to the diagonally opposite one. Taking 
both of these paths better maintains the symme- 
try of the ground state potential and therefore 
provides improved overlap with the ground state. 
We also consider 2x3 boxes which are the aver- 
ages of the possible paths connecting two opposite 



Figure 3. Projection of the Y-shape path on the 
xy-plane. 



corners using 1x1 and 1x2 boxes. We will create 
further link paths in the future for use in bigger 
loops as necessary. Hence a diagonal staple is, in 
fact, an average of several "squared path" staples 
connecting the same end points. 

The quark co-ordinates considered for Y-shape 
paths are summarized in Tabled We note that 
the origin (0, 0) is not at the centre of these co- 
ordinates. Rather the coordinates are selected to 
place the quarks at approximately equal distances 
from each other. 

3. GLUON FIELD CORRELATION 

In this investigation we characterize the gluon 
field by the action density S{y, t) observed at spa- 
tial coordinate y at Euclidean time t measured 
relative to the origin of the three-quark Wilson 
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Figure 4. The diagonal path (sohd hne) is taken 
as the average of the two possible ones going 
around the box enclosing the end points. 



(x, y) Coordinates Separation 
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Cart. 


Diagonal 


(0,1) 


(1, 


-1) 


(-1,-1) 
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2.24 


(0,2) 


(2, 


-1) 


(-2,-1) 
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3.61 


(0,3) 


(2, 


-1) 


(-2,-1) 
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4.47 


(0,3) 


(3, 


-2) 


(-3,-2) 
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5.83 


(0,4) 


(4, 


-3) 


(-4,-3) 
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8.06 


(0,5) 


(5, 


-4) 


(-5,-4) 


10 


10.3 



Table 1 



(x, y) coordinates for the three quarks consid- 
ered for Y-shape paths. The separation of the 
two quarks along a Cartesian (Cart.) direction 
are compared to the separation of those from the 
third quark. 



loop. We calculate the action density using the 
highly- improved C(a'*) three-loop improved lat- 
tice field-strength tensor JH] on four-sweep APE- 
smeared gauge links. 

Defining the quark positions as fi, r2 and ra 
relative to the origin of the three-quark Wilson 
loop, and denoting the Euclidean time extent of 
the loop by T, we evaluate the following correla- 
tion function 

C(y;ri,f2,f3;T) = 

{W:,Q{r^,f2.r^-T)S{y,Tl2)) 
(W^3Q(ri,r2,f3;T))(5(y,r/2)) ' 

where (• ■ ■ ) denotes averaging over configurations 
and lattice symmetries as described below. This 
formula correlates the quark positions via the 
three-quark Wilson loop with the gauge-field ac- 
tion in a gauge invariant manner. For fixed quark 
positions and Euclidean time, C is a scalar field in 



three dimensions. For values of y well away from 
the quark positions r^, there are no correlations 
and C ^ 1. 

This measure has the advantage of being pos- 
itive definite, eliminating any sign ambiguity on 
whether vacuum field fiuctuations are enhanced 
or suppressed in the presence of static quarks. 
We find that C is generally less than 1, signaling 
the expulsion of vacuum fluctuations from the in- 
terior of heavy-quark hadrons. 

4. STATISTICS 

For this work we consider 200 quenched QCD 
gauge-fields using the C'(a^)-mean-field improved 
Luscher-Weisz plaquette plus rectangle gauge ac- 
tion 22 on 12^^x24 lattices at /3 = 4.60, providing 
a lattice spacing of a = 0.123(2) fm as set by the 
string tension. 

To improve the statistics of the simulation we 
use various symmetries of the lattice. First, we 
make use of translational invariance by comput- 
ing the correlation on every node of the lattice, 
averaging the results over the four-volume. 

To further improve the statistics, we use re- 
flection symmetries. Through reflection on the 
plane a; = we can double the number of T and 
Y-shaped Wilson loops. By using reflections on 
both the plane x = Q and y = we can quadruple 
the number of L-shaped loops. 

We finally use 90° rotational symmetry about 
the a;-axis to double the number of Wilson loops. 
This means we are using both the xy and xz 
planes as the planes containing the quarks. 

For this exploratory study we present images 
for r = 2 and note that qualitatively similar re- 
sults are observed for T = 4. We plan to examine 
this and other similarly conservative alternatives 
on our anticipated larger lattice volumes provid- 
ing further statistics improvement. 

5. SIMULATION RESULTS 

Figure plots the correlation function of 
Eq. C{y), for an ortho-slice in the plane of 
the three-quark system. The colours of the ortho- 
slice, denoting the scalar values of C{y), are ren- 
dered as a surface plot below. Away from the 



Figure 5. Ortho-slice and associated surface plot 
of the correlation function of Eq. C(y), ob- 
served for positions y in the plane of the three- 
quark system. 



quark positions C{y) — > 1, but within the three- 
quark system C{y) ~ 0.9 indicating the sup- 
pression of QCD vacuum field fluctuations from 
the region inside heavy-quark hadrons. Here the 
quarks are approximately 0.75 fm from the cen- 
tre of the distribution and are about 1.25 fm from 
each other. 

Figure El renders the spatial points of C{y) 
where the gluon action density is largely sup- 
pressed. Tube-like structures are revealed in Y- 
shape as opposed to A shape in accord with ex- 
pectations from precision static quark-potential 
analyses. 

We note that the centre of the flux tube is 
equidistant from the quark positions and not cen- 
tered over the origin of the Wilson loop. This 
gives us confidence that sensitivity to the source 
is minimal. However it is important to examine 
this issue in greater detail and this will be the 
subject of a forthcoming publication. 
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